Bistable hysteresis and resistance switching in hydrogen-gold junctions 
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Current-voltage characteristics of H2 — Au molecular junctions exhibit intriguing steps around 
a characteristic voltage of V 3 « 40 mV. Surprisingly, we find that a hysteresis is connected to 
these steps with a typical time scale > 10 ms. This time constant scales linearly with the power 
dissipated in the junction beyond an off-set power P s = IV S . We propose that the hysteresis 
is related to vibrational heating of both the molecule in the junction and a set of surrounding 
hydrogen molecules. Remarkably, by stretching the junction the hysteresis' characteristic time 
becomes > days. We demonstrate that reliable switchable devices can be built from such junctions. 



The central idea behind molecular electronics is that 
novel devices can be created using intelligently de- 
signed functional molecules. The simplest example is 
a molecular junction that exhibits conductance switch- 
ing. Chemists have synthesized various types of bistable 
molecules, many of which have subsequently been applied 
in devices 1 ^ 2 -' 3 -'^. Surprisingly, however, conductance 
switching has also been observed in junctions that incor- 
porate passive molecules only^^^. Understanding the 
fascinating physics behind the latter phenomenon may 
lead to a novel design philosophy of switchable molecu- 
lar devices. Here, we present detailed transport studies 
on H2 — Au molecular junctions. We demonstrate that 
if I(V)-curves are taken by fast data acquisition, hys- 
teretic behavior is observed. Moreover, we can manipu- 
late our junctions such that the hysteresis' characteristic 
time scale grows from ~10 ms to >days. The resulting 
H2 — Au junctions are arguably the smallest switchable 
devices around. The physics behind these fascinating 
phenomena forms the main subject of this Letter. 

Molecular junctions incorporating hydrogen molecules 
have a number of remarkable propertie s 10 ' 11 ' 12 ' 13 ' 14 . 
Perhaps most striking is that differential conductance 
(dl/dV) curves exhibit sharp features at a well-defined 
voltage K' 11 ' 12 ' 13 ' 14 - In FiglU we show two typical mea- 
surements for a H2 — Au contact. In one case, two sharp 
peaks are observed (FigQJ,), whereas in the other, two 
dips are seen (FigQJ)). Both curves are symmetric with 
respect to voltage. Remarkably, the dip in FigHt) is so 
deep that the differential conductance around V s = ±43 
mV becomes negative. Recently, Thijssen and Halbritter 
et al. measured curves similar to Fig. [T£i) and b) and 
put forward a two- level mode l 11 ' 12 . The two levels (1, 
2), most likely representing different geometries of the 
H2 — Au contact, give rise to two specific conductance 
values (G±, G'2) 1 ^. Hence, if transitions between level 
1 and level 2 are induced at V = V s , the conductance 
changes discontinuously from G\ to an average of G\ and 
G2 (the weighing being dependent on the degeneracy of 
the levels^ 2 -). This naturally results in a dip (if G\ > G2) 
or peak (if G\ < G2) in dl/dV, or analogously, in a step 
down or up in an I(V)-curve. The latter is illustrated by 
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FIG. 1: a),b) Examples showing peaks (a) or dips (b) in 
dl/dV around w 43 mV, measured on a H2 — Au contact 
(at 5 K). Note that (b) features a negative differential con- 
ductance, c) I(V)-curves measured with varying sweep rates. 
The lower curve, measured in 10 s, corresponds to curves such 
as a). However, for faster voltage sweeps, a hysteresis shows 
up around V s . Sweeping times are 10, 1, 0.5, 0.33 and 0.1 
s. For clarity, the traces are shifted vertically, d) I(V)-curve 
showing both hysteresis and two level fluctuations around V 3 
(voltage swept in Is). Inset: double well potential for the two- 
level model. We associate the fast switching between the two 
wells to time scale r\, and the slower changes in the potential 
landscape to ti. 



the lowest I(V)-characteristic presented in Fig. Q}:). In 
both models proposed, a double potential well is assumed 
to describe the two- level system (schematized in FigQJl) , 
inset). However, there is disagreement on the question 
how a transition between the levels occurs. Halbritter et 
al. relate V s to the energy difference between the two lev- 
els. Furthermore, they explain the negative differential 
conductance by proposing that level 2 is highly degener- 
ate. In contrast, Thijssen et al. propose that a transition 
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is only possible if a molecular vibration is excited with 
high enough energy fkJi to assist a jump over the bar- 
rier. Hence, they conclude that V s — Hu>i/e, a statement 
for which they provide significant experimental evidence. 
Regardless of the exact model, however, the transition 
between states 1 and 2 is anticipated to be fast around 
V s , with time scales < 1 ns, determined by tunnelling or 
electron-phonon interaction. The results we present here, 
strongly contradict this expectation. In FigQJ), we show 
I(V)-curves taken with various sweeping times. When 
slowly ramping V, one obtains the lower curve discussed 
above. However, when sweeping faster than 1 s per curve, 
a pronounced hysteresis appears around V = V S . In other 
words, there is a macroscopic time scale connected to the 
system's return from the high conductance state to the 
low one. Moreover, for even faster voltage sweeps, > 10 
Hz, the system does not have enough time to relax back 
to state 1; it keeps the 'high' conductance value for all V. 
Not only does this not agree with the two-level models, it 
also rules out the possibility that the steps in the I(V)'s 
are due to resonant tunnelling. Hence, this intriguing 
hysteresis forms the central item of this study. 
To obtain stable molecular junctions, we employ mechan- 
ically controllable break junctions^. Basically, a gold 
wire is broken in cryogenic vacuum (at «5 K), and the 
electrode tips are sharpened by a " training" procedure^. 
Subsequently, a small volume of Hi (99.999 % pure) is 
introduced and the junction is opened with «5 pm/si^. 
Meanwhile, I(V)-curves are recorded by sweeping V be- 
tween —90 and 90 mV at 1 Hz. After several minutes, 
nonlinear I(V)-curves are observed with pronounced con- 
ductance step o 11 ! 12 . During this procedure, we do not 
observe a jump out of contact. In total, we have used 6 
break junctions on which 56 different Hi — Au contacts 
were studied. To obtain geometrically different Hi — Au 
junctions for one single device, we use the following pro- 
cedure. First, the Hi — Au junction is broken, result- 
ing in a jump out of contact. Then, the electrodes are 
pushed back into contact to conductance values > 10Gn 
(where Go = 2e 2 /h). Finally, the gentle pulling process is 
started again until a new Hi — Au contact forms. To con- 
nect to previous work, we have first investigated the I(V) 
curves of all these junctions at low voltage scan rates. 
In Fig. we present a statistical analysis of the on- 
set voltage V s . Interestingly, V s is positioned around 40 
meV, rather independent of the conductance of the con- 
tact. Indeed, the zero bias conductances show a broad 
variation between 0.003 and 0.3 Go- Since we did not ob- 
serve a jump out of contact while creating the junctions, 
this may point to chain formation as reported in RefJ^. 
These authors find conductance values in the same range 
as observed here. In Fig. [2b), we present a histogram 
of all V s values. The average position is V s — 40 ± 3 
mV. This is in good agreement with Thijssen et al, who 
found V s = 42 ±2 mV— . As explained above, they relate 
Vg to a Hi phonon, which induces transitions between 
levels 1 and 2. Experimentally, they provide two pieces 
of evidence. First, they show that V s shifts down by a 



factor w y/2 when Hi is replaced by Di (deuterium). 
Second, they show that upon pulling a junction, a step 
down in a dl/dV-curve can evolve into a peak. Since a 
step in dl/dV is related to the excitation of a phono n 20 ' 21 , 
this confirms their assertion. We have performed similar 
pulling experiments. In some cases, peaks in dl/dV grew 
larger. In other cases, they would virtually disappear, 
leaving a step as the dominant feature. This confirms 
the result by Thijssen. 
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FIG. 2: a) Position of the conductance peak/dip for 56 dif- 
ferent contacts as a function of its zero bias conductance, b) 
Histogram of all data points. On average, V a — 40 ± 3 m V. 

As we increase the scan velocity of our I(V)-scans, a 
hysteresis develops around V s (Fig[Tb). Generally, we ob- 
serve a smooth transition from one conductance state to 
the other. However, this is not always the case. The 
I(V)-scan in Fig[TJi) for example, shows pronounced two 
levels fluctuations (TLFs) in combination with hystere- 
sis. A comparison of FigsQJi) and d) leads to a number of 
important observations. First, the hysteresis is not an ar- 
tifact of our electronics: TLFs much faster than the time 
scale of the hysteresis are easily observed. Second, the 
occurrence of TLFs confirms a two-level model. Halbrit- 
ter et al. could not resolve the actual TLFs. Also in 
our case, most curves show an apparently smooth step 
around V s . This is due to the TLF frequency exceeding 
the bandwidth of the electronics (« 50 kHz). Third, we 
see that the TLF-amplitude around Vs is approximately 
equal to the difference between the 'low' and 'high' con- 
ductance states far away from V s . This indicates that for 
V ^> V s , the system is in state 2 most of the time. Hence, 
state 2 is either relatively stable or highly degenerate! 2 ". 
Finally, we deduce from Fig. [TJi) that despite contin- 
uous communication between both levels (through the 
TLFs), the hysteresis is still present. This suggests that 
the exact shape of the double potential well (at V » V s ) 
differs for the upward and downward voltage scans. We 
deduce that there are two time scales in our system: a 
time Ti corresponding to fast switching between states 
(TLFs) and ti > n connected to changes in the effective 
potential landscape, resulting in hysteresis. 

To gain more insight in ti, we have recorded response 
times using well-defined voltage pulses. An example is 
given in Fig[3k.). First, V is set to a value Vhigh > V s , i.e. 
beyond the conductance step. The system's conductance 
is 'high'. After 1 s, long enough to reach equilibrium, the 
voltage is instantly dropped to a value Vi ow < V s (in this 
case 20 mV). Interestingly, the current responds to this 
decrease within 20 us, consistent with our bandwidth. 
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However, the conductance I/V does not. It is still 'high' 
after 20 /is. Only after 10 ms, it starts to drop to a lower 
value. We define t 2 as the time in which the conductance 
drops halfway (dotted line in FigJSJa). Remarkably, we 
find larger t 2 for higher Vhigh (see Fig|3p), which implies 
that T2 depends on the history of the junction. We have 
systematically investigated t 2 vs. Vhigh for 4 different 
contacts, as depicted in Fig|3b)2£. Clearly, r 2 increases 
more than linearly with Vhigh, with an onset around 40 
meV. In fact, if we plot r 2 as a function of I ■ {Vhigh — V s ) 
we find a linear relation for all 4 samples (see Fig. [3Ji). 
Hence, there is a direct connection between t 2 and the 
power IVhigh that was being dissipated before the voltage 
step, corrected for an onset value P s = IV S . From this 
interesting result, we deduce that heating effects (multi- 
ple phonon excitations) are at the basis of the hysteresis. 
The fact that V s = 40 mV serves as a threshold in all 
cases points to the pivotal role of a hydrogen phonon in 
this process^. A second experiment that demonstrates 
this connection is shown in Fig. (3^)- This measurement 
is similar to that in Fig. [SJd). However, now Vhigh is 
kept constant at 80 mV, whereas V ow is systematically 
varied. We find that r 2 increases dramatically as Vi ow 
approaches V s . 
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FIG. 3: a) Pulse response. At t=0, the voltage is instantly 
dropped from a high voltage (Vhigh > V s ), to a low voltage 
(Viou, < V s ). Here, Vhigh (Vi ow ) is 50 mV (20 mV). Right 
after the voltage drop, the conductance is still in the high 
conductance state, b) Response for different initial voltages. 
Vhigh is systematically varied between 50 mV (first curve) and 
90 mV (last curve), while Vi ow is 20 mV. The decay time 
T2 is defined as the time where the conductance has dropped 
halfway (16 ms for Vhigh =50 mV). c) Decay times as a func- 
tion of Vhigh, for 4 different contacts (Vi ow = 20 mV.) The 
contacts have a conductance of, from lower to upper curve, 
0.015, 0.044, 0.048 and 0.058 Go (at 50 mV). d) Decay times 
as a function of dissipated power during the voltage pulse, cor- 
rected for an onset value P s = IV 3 (same contacts as in c). 
e) Decay times as a function ofVi ow . Vhigh was set to 80 mV 
and for this contact V s = 39 mV. 

Remarkably, when slowly stretching the contact, a spe- 
cial type of hysteresis appears for the majority of the 
junctions. In general, a contact is stretched in steps of 



10 pm while for each position an I(V) curve is measured 
(in 1 min). Figf^t) shows 4 I(V)'s with 50 pm intervals. 
The upper curve is measured right after the H 2 — Au con- 
tact is formed. With further stretching a hysteresis ap- 
pears, which is most pronounced after 2 A (lower curve). 
Interestingly, this type of hysteresis does not show any 
time dependence, contrasting the data above. Therefore, 
we can use it as a stable memory element. For this, 
we first apply an offset voltage V = V s . Subsequently, 
we switch conductance by applying positive (low G — > 
high G) or negative voltage pulses (v. v.). This is demon- 
strated in FigfJ}). The switch is remarkably stable in 
time. Even after 14 hours, no degradation is observed 
(see Fig. Ht). Furthermore, no deterioration is seen after 
multiple switching events (see Fig. Hli). We have per- 
formed over 5000 switching events, and saw no changes to 
the sample. Hence, we have created a switchable molec- 
ular device from the smallest molecule nature provides. 
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FIG. 4: a) I(V) sweeps while pulling the electrodes (voltage is 
swept in 1 min). After stretching 2±0.3j4, a time independent 
hysteresis appears. For clarity, the curves are shifted verti- 
cally. Insets: two other examples of contacts showing time 
independent hysteresis, b) Hydrogen switch. Upper graph: 
Voltage used to control the switch. On top of an offset voltage 
(here 34 mV), the voltage is pulsed with +20 mV and -20 mV 
(duration is 20 ms) . Lower graph: Pulse response. After each 
pulse, the system switches, c) Stability of the hydrogen switch. 
After 10 hours, the voltage is pulsed by +20 mV (20 ms) to 
switch to Gi- After 14 hours, the measurement is stopped, d) 
Multiple switching events. The pulse height is ±4 mV, on top 
of a dc voltage of 34 mV (pulse duration is 20 ms). 

To discuss our observations, we first return to Figj3] 
Here we concluded that heating, i.e. electron-phonon ex- 
citation, plays a pivotal role in the hysteresis. Further- 
more, the molecular phonons are crucial in this process. 
We never observed heating effects on bare gold electrodes, 
which is consistent with the relatively low dissipation in 
our junctions (cf. Ref.— ). Recently, Galperin et al. and 
Chen theoretically studied the temperature of molecu- 
lar junctions under a bias voltag o 23 ' 24 . It is not trivial, 
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however, to define a temperature for a single molecule. 
Galperin et al approach this problem as follows. First, 
they assume that an external phonon bath is weakly cou- 
pled to the active mode of the molecule in the junction. 
Then, they calculate the ensuing heat flux between mode 
and bath. The bath temperature at which the heat flux 
vanishes is finally assigned to the molecular junction. For 
^-junctions, a temperature increase ~ 100 K was pre- 
dicted, exactly when the voltage reaches the phonon en- 
ergy eV = 7ta. i 23 ' 24 . From this, one might naively expect 
a hysteresis between upward and downward I(V)-scans. 
Going up in V, the molecule is 'cold' and in state 1; com- 
ing down, the molecule is 'hot', i.e. in a combination of 
states 1 and 2. However, it is highly unlikely that a sin- 
gle molecule will remain excited on time scales as large 
as T2 > 10 ms. For the most likely relaxation mecha- 
nism, electron phonon coupling, typical time scales are 
< 1 ns, depending on the vibration mode 25 . We pro- 
pose that there is a larger phonon bath that preserves 
the energy for a relatively long time (cf. Galperin). A 
good candidate is formed by the large set of H 2 molecules 
that likely surround the actual molecule in the junction. 
At V > V s , these molecules will be excited by efficient 
phonon-phonon coupling from the junction outwards. Af- 
ter V is dropped below V s , it will take a relatively long 
time before this bath (and the molecule in the junction) 
is at the base temperature again. It may well be that 
heating of surrounding H2 molecules actually induces a 
phase transition, although this is not essential to explain 
Figs. Q][3] In fact, we have only observed current steps at 
experimental temperatures below « 20 K. A phase tran- 
sition has been proposed by Gupta and Temirov et al. to 
explain non-linear IV curve o 13 i 14 . Such a transition could 
indeed explain the time-independent hysteresis observed 



in Fig. [4] Upon a phase transition, the hydrogen phonon 
energy will most likely change. If it assumes a lower value 
u>2 than before (wi), the hydrogen molecule (s) will still 
be excited if the voltage drops below V s , provided that 
LU2 < eV < ll>\. This situation may be quasi-permanent, 
since the system is continuously excited. Only when V 
is decreased further, a transition back is expected. Our 
work may thus represent the first experimental case in 
which molecular heating is directly distinguished in a sin- 
gle junctio n 23 i 24 . 

In summary, we demonstrate that the presence of H2 
in a Au junction leads to two-level fluctuations, steps 
in I(V)-curves and hysteresis. Moreover, by carefully 
stretching the junctions, a reversible and stable molec- 
ular switch can be created. Interestingly similar fea- 
tures have been reported for a large variety of molec- 
ular junctions, with molecules ranging from larger or- 
ganic molecules to CO^ ' 11 ' 26 ' 27 . Since such switching 
effects are general phenomena in molecular electronics, 
there may be one basic mechanism behind them. In 
parallel, our work stresses the requirement to work in 
a hydrogen-free environment if one focuses on a differ- 
ent molecular system. Finally, we note that the steps in 
I(V)-curves (explained by a two-level model) may easily 
be confused with resonant tunnelling. To distinguish the 
two, fast I(V)-scanning is an efficient tool, since resonant 
tunnelling should not lead to hysteresis. 
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